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The mesoscale rheology of a lamellar phase of a block copolymer is modeled as a 
structured fluid of uniaxial symmetry. The model predicts a viscoelastic response 
that depends on the angle between the the local lamellar planes and velocity gradi- 
ents. We focus on the stability under oscillatory shear of a two layer configuration 
comprising a parallel and a perpendicularly oriented domain, so that the two layers 
have a different viscoelastic modulus G*{u}). A long wave, low Reynolds number 
expansion is introduced to analytically obtain the region of stability. When the re- 
sponse of the two layers is purely viscous, we recover earlier results according to which 
the interface is unstable for non zero Reynolds number flows when the thinner layer 
is more viscous. On the other hand, when viscoelasticity is included, we find that 
the interface can become unstable even for zero Reynolds number. The interfacial 
instability is argued to dynamically favor perpendicular relative to parallel orienta- 
tion, and hence we suggest that the perpendicular orientation would be selected in a 
multi domain configuration in the range of frequency oj in which viscoelastic contrast 
among orientations is appreciable. 

I. INTRODUCTION 

A stability analysis of two superposed fluid layers to oscillatory shear is given when the 
frequency dependent viscoelastic modulus of each layer is different. At vanishing Reynolds 
number, we find that an initially planar interface can be rendered linearly unstable to long 
wavelength perturbations by viscoelastic contrast. Our analysis is motivated by rheology and 
alignment studies of structured phases, mostly in soft and biological matter, in which local 
viscoelastic response couples to a broken symmetry variable of the phase (e.g., orientation). 
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The specific case considered is related to alignment studies of lamellar phases in bulk block 
copolymer melts. There by reason of symmetry, the complex viscoelastic modulus G*{uj) 
depends on the relative orientation between local lamellar normal and velocity gradient. 
We argue that given the dependence of G*{u) on orientation, the hydrodynamic instability 
discussed here would dynamically favor domains of perpendicular orientation over parallel 
in a bulk, multi domain configuration. 

Block copolymers are currently being investigated for a wide variety of applications in 
nanotechnology^"^ due to their ability to self-assemble into ordered nanophases of different 
symmetry (lamellar, cylindrical, spherical, close packed spherical, or bi-continuous phases 
such as gyroid).'"'' Equilibrium morphologies and characteristic length scales can be easily 
controlled by chemical manipulation of the polymer blocks. In general, however, quenched 
block copolymer samples are macroscopically disordered in that they contain a large number 
of domains of different orientation, all degenerate by symmetry. Given that disordered or 
defected samples are not generally suitable for applications^^"^^, a number of strategies 
have been put forward to accelerate annealing of defects, or to post process polycrystalline 
samples in order to increase the characteristic size of ordered domains. Applied shears have 
been particularly effective in aligning thin films of in-plane cylinder sphere forming^^, and 
lamellar^ phases. The same strategy has been used on bulk samples^^"^^, although there is no 
explanation yet as to the selected orientation relative to the shear as a function of material 
parameters, and amplitude and frequency of the shear^^. This lack of understanding of the 
non equilibrium processes involved in the alignment of bulk samples has led to alternative 
strategies including, for example, solvent annealing and applied electric fields^". 

Our current understanding of shear alignment in bulk samples is well summarized by Wu 
et al.^^. They conducted a systematic analysis in styrene-isoprene (SI) multiblocks (from 
diblocks to undecablocks). Results for diblocks confirm earlier experimental findings that 
show a sequence of transitions from parallel orientation at low frequencies, to perpendicular, 
and again to parallel as the frequency is increased. The sequence of transitions is quali- 
tatively different in poly- (ethylene-propylene) - poly-(ethylethylene) (PEP-PEE) diblocks, 
presumably because a marked viscoelastic contrast between styrene and isoprene blocks 
which is largely absent in PEP-PEE. This issue was already emphasized by Fredrickson^^ 
who noted that styrene is largely unentangled, whereas isoprene is entangled, and hence 
a large contrast in relaxation times of the blocks is anticipated. Further, the experiments 
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of Wu et al.^^ show that this double transition is observed in diblocks and triblocks, but 
surprisingly not in higher order multiblocks, in which entanglement effects among longer 
chains would be expected to be more important. 

A further complicating factor in the elucidation of the orientation selected under shear is 
the non terminal viscoelastic response of the block copolymer melt in the range of frequencies 
in which orientation switching is observed. ^^'^^"^^ (terminal behavior has an elastic modulus 
that vanishes with frequency as G' ~ uj"^ whereas the loss modulus reflects Newtonian viscous 
dissipation G" ~ w). In fact, non terminal rheology is observed even at the lowest frequencies 
at which G*{uj) has been measured, a range well below that which can be attributed to 
individual chain dynamics. Interestingly, a dependence of the non terminal viscoelasticity 
on lamellar orientation has also been found in largely oriented samples. ^^'^^"^"^ Except for 
diblocks, all multiblocks in Ref. 19 show terminal behavior as w — ?■ when aligned along 
the perpendicular orientation, whereas none of them are terminal in the parallel orientation. 
This observation, together with the prominent entanglement effects described above, calls 
into question existing theoretical analyses that either neglect hydrodynamics or, if they 
don't, they assume Newtonian flows. Viscoelastic response appears to be closely correlated 
with orientation, and is certainly not negligible even at the lowest frequencies probed. 

We propose a new viscoelastic model of block copolymer melts that explicitly takes into 
account the fact that the melt in a lamellar phase is a uniaxial (not isotropic) fluid. As such, 
local viscoelastic response involves the coupling between local orientation of the lamellae 
and local strain and velocity gradients. For an incompressible system, such a symmetry 
requires three viscosity coefficients and three elastic moduli^®. We extend here the earlier 
analysis of a purely viscous uniaxial fluid^^, to incorporate viscoelasticity. As shown below, 
a uniaxial viscoelastic medium readily allows for differential rheology in, say, parallel and 
perpendicular orientations, including the observed terminal behavior of the perpendicular 
orientation versus viscoelastic response of other orientations. 

For simplicity, our analysis does not address a macroscopically disordered sample, but 
rather a configuration comprising only two perfectly oriented domains and the boundary 
separating them. Of the three possible orientations of an ordered lamellar domain relative 
to the imposed shear we focus here on parallel and perpendicular only. The third indepen- 
dent component (the so called transverse, in which the lamellar normal is parallel to the 
velocity direction), is less stable as it is being compressed by the shear^*^. Both parallel and 



4 



perpendicularly oriented domains are marginal with respect to the flow, and as such have 
been the subject of most orientation selection studies to date. 

Section II describes the model equations and the fluid configuration under consideration. 
We introduce a general linear viscoelastic constitutive law for a incompressible system with 
uniaxial symmetry, and derive the reference flow under oscillatory shear and a planar bound- 
ary separating two layers of different orientation. Section III presents our results concerning 
the linear instability of the base state by considering small perturbations of large wavelength, 
also in the limit of small Reynolds number. The implications of our findings on orientation 
selection in a lamellar block copolymer are discussed in Sec. IV. 



Although the description below can be generalized to the case of a multi block (see, e.g., 
Ref. 6), we start by introducing the two fluid model of a diblock copolymer for monomers 
of type A and B. At sufficiently low frequencies so that the individual polymer chains have 
relaxed to their equilibrium conformation, the state of the copolymer can be described by a 
local order parameter ip ~ Pa—Pb, where pA and ps are the number fractions of monomers A 
and B respectively. The evolution of the order parameter has a relaxational component that 
is driven by free energy reduction^^ If fluid flow is allowed, there is a reversible contribution 
to the order parameter equation that incorporates advection by the local flow. If Vi is 
the i — th component of the local fluid velocity, the advection term for an incompressible 
fluid is Vidiip, where sum over repeated indices is assumed. The flow velocity Vi satisfies 
the momentum conservation equation with the stress tensor being the momentum density 
current. We now introduce a general linear constitutive law for the stress tensor of the form. 



The independent contributions to a fourth rank tensor compatible with incompressible con- 
dition {6ijdiVj = 0) and uniaxial symmetry (gj is the i—th component of the local unit normal 
to the lamellar planes, and assume that the constitutive law is invariant under qi — — g^) 
can be written as 



II. VISCOELASTIC MESOSCOPIC MODEL 




(1) 



Gijkiif) =Gi{t)qiq,jqkqi + Gi{t){5ik5ji + 6ii6jk) 

+ G^eit) qi{qkSji + qi6kj) + qjiqkSn + qiSki) , 



(2) 
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where Sij is the Kronecker's delta. There are, in this case, three independent modu- 
lus components Gi,G4 and G56 and, as a consequence, the viscoelastic response depends 
on the local lamellar orientation. For example, for perfect parallel lamellae (q || Vv) 
G\\{t) = Giit) + G^Qit), and for perfect perpendicular lamellae (q || V x v) G±{t) = G^it). 
Therefore, any viscoelastic contrast between these two orientations is contained in the com- 
ponent G56(t). Note, in particular, that if the response of the material in the parallel 
configuration is terminal, this has to be the case for a perpendicular orientation as well. 
The reverse, however, is not true, fact that is consistent with experiments^^. 

We focus here on extending the results of Ref. 27 to include uniaxial viscoelasticity. 
We consider two semi infinite layers of a lamellar block copolymer, one on top of the other 
(Fig. 1), confined between two infinite and parallel planes. The phase at the top has thickness 
H, and the lower phase h. The upper plane oscillates with a velocity U of frequency u and 
amplitude U while the lower plane is at rest. The densities of two layers, p, are constant 
and equal since they are the same copolymer, differing only on lamellar orientation. For 
simplicity, we neglect any variation in the vorticity direction (V x U), ot y in Fig. 1. 
Hence we focus on a two dimensional system in which the fluid velocity has components 
= (u*-"-*, 0, u;'-"'') with a = 1,2 denoting fluid 1 (lower) or 2 (upper). 

When the two layers in this configuration contain parallel and perpendicular phases, 
the state of the system is unaffected by the flow since v ■ Vip = for both orientations. 
Furthermore, if as is customary, both phases are assumed to be isotropic fluids, then the 
configuration shown is stable, the order parameter is stationary, and the flow is a simple 
linear shear. If, on the other hand, Eq. (2) is assumed, then the two layers have a different 
viscoelastic modulus and the configuration can be linearly unstable. This instability is the 
subject of our study below. 

The response of stratified fluids to shear flows has been studied in detail. Following 
pioneering work on the interfacial stability of two superposed, incompressible and Newtonian 
fluids under steady and oscillatory shears by Yih'^°''^^, there have been a series of additional 
studies involving viscous fluids under steady^^''^'^ and oscillatory shears^'^, or viscoelastic 
fluids under steady shear^^"'^^. To our knowledge, however, there has been no stability 
analysis of stratified viscoelastic fluids subjected to oscillatory shears. In the case of purely 
viscous fluids, a stratified configuration can be unstable because of a phase lag between 
the velocity fields on either side of the interface. For this to be possible in a Newtonian 
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FIG. 1: Schematic configuration studied. As the top and bottom fluids represent lamellar phases 
of parallel or perpendicular orientation, the complex modulus G*^°'^{uj) of the two is different. The 
distortion of interface away from planarity is given by Zg = Zs{x,t). 



fluid, the Reynolds number must be finite. This case is unhkely to be of relevance to block 
copolymers as typical Reynolds numbers are extremely small. In the viscoelastic case, on 
the other hand, a differential elastic response between the two phases, even at zero Reynolds 
number, provides for a mechanism for instability as detailed below. 
In each fluid layer, the velocity field satisfies 



where the stress tensor is 



(3) 



(4) 



where we have separated P*^"-* , the hydrostatic pressure, from the stress a^"^ . Incompressible 



flow is assumed in both layers 



dJ''^ = 0. 



We now consider a different viscoelastic modulus in each layer and write 



a. 



(a) 



dt' G^''\t - 1') \d.jv'f\t') + div^;'\t') 



(5) 



(6) 



In order to make the problem as general as possible, we do not specify any functional form 
of the moduli of the two viscoelastic fluids, although we must point out that our study is 
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confined to linear viscoelasticity. We next define the complex modulus^^ 

POO 

G*{uj) =iu dt G(t)e-*^* = G'{u) + iG"{uj), (7) 



where is the storage modulus, and G"{u) the loss modulus. 

The interface separating the two layers is Zg = Zs{x, y, t). Since fluid particles move with 
the interface (no mass transport across the interface is allowed) we have, 

5t^s + w^°^|.=..9,2s = w;(°)|,=,,. (8) 

The set of governing equations (3) - (5) is to be supplemented by boundary conditions. 
We use the no-slip and no-penetration boundary condition on the planes. For the upper 
fluid, 

^ {z = H,t) = {U cos{ujt) ,0,0), (9) 

and, for the lower fluid, 

vl^\z = -h,t) = 0. (10) 
At the interface {z = Zg) the velocity is continuous 

V?\z = Zs)=vi'\z = Zs). (11) 

The stress tangential to the interface is continuous whereas the stress normal to the interface 
is discontinuous and balanced with the interfacial tension T 

7ilfnjti\,=,^ = 7iifnjti\,=,^, (12) 

7rlfnjni\,=,^ - 7i^fnjni\,=,^ = -T{dl + dl)zs, (13) 

where Ui and tj are the unit vectors normal and tangential to the interface. Therefore, 
Eqs. (3) - (13) completely describe the dynamics of two superposed viscoelastic fluids under 
oscillatory shears. 

A. Reference state 

In the base reference state the interface is planar {zs = 0), and the flow is stratified and 
parallel to the planes, with a gradient along z-axis, and oscillating with the driving frequency 
u. Hence v^""* = u^^\z,t)x and we assume 

u'^^\z,t) = u^^\z)e-"^' + ^(z)e*"*. (14) 
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When v["'' is substituted into the momentum conservation equation, its ^-component be- 



'b 



comes diP^°'^ = 0, leading a constant hydrostatic pressure ■ On the other hand, the 
equation for the x-component involving the mode e"**^* is 



f dt' G^°'\t-t')e-'''''dlu^^\ (15) 

J — oo 



(a) r (a) 

-= \u_ ' 



For the e"^* mode, m+ satisfies the complex conjugate of Eq. (15); therefore, u 
and M^""* is real. With the definition of the complex modulus Eq. (7), Eq. (15) becomes 

Note that, since G{t) is real, G*{-uj) = \g*{uj)Y = G'{uj) - iG"{uj). 

The boundary conditions for the Fourier e~*"* mode are m^'' {z = H) = U/2 and u'-}\z = 
—h) = 0. At the planar interface rij = (0,0, —1), ti = (1,0,0) so that continuity of velocity 
is u^^\z = 0) = u^}\z = 0). In addition, the continuity of tangential stresses at the interface 
leads to a pressure that is continuous and constant across the interface Pj^"^ = P^. The 
balance condition of the normal stresses at the interface provides the following restriction 
to fiuid motion in the layers 

G<^\-u)d,u^^\z = 0) = G<^\-uj)d,u-l\z = 0). (17) 

Since the interface is not distorted in the base state, the effect of the surface tension is 
absent. 

We now solve the differential equation (16) by assuming 

u^^\z) = A^^^^"''^ + fi(")e-'='"'^ (18) 

Upon substitution into Eq. (16), and using all four boundary conditions, we find 

[z) = ^ sinh k^^\h + z), (19) 
2K 



(2)/ \ ^ 



sinh k^^^h cosh k'^'^^z + ^ — ^} , /n^ cosh k^^^h sinh k^'^'^ 

G<^){~uj)k^^> 



'z 



(20) 



where 



and 

K = sinh k^^^h cosh k^^'^H + -—^^ ^ cosh A; ^/i sinh A; if. 

G*i^>{-u)ky^) 



(22) 
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Note that the reference state exphcitly depends on the complex moduli at the driving fre- 
quency oj. In the next section we study the linear stability of this base flow against a small 
perturbation of the interface. 

III. LINEAR STABILITY ANALYSIS 

We introduce small perturbations away from the base state v = v;, + ^v^"'', with ^v^"-* = 

Since our system is effectively two dimensional, 
we introduce the stream function 0*^"^ such that 5u^°'^ = dz4>^°'^ and 5w^°'^ = —dx4>^°'^. We 
decompose the perturbations into normal modes in the x direction 

0(a) = ^p{") = 5pH(2^ t)e^9.^ (23) 

and 

Zs = UtY'''^ (24) 

In general, (^'^"■'(z, t), Zs{t) and SP^"\z,t) are complex and depend on q^. The resulting 
differential equations, after linearizing in the amplitude of the perturbations, have time 
periodic coefficients in 27t/u, and hence we use Floquet theory. We therefore write 

= 0(")(2,t)e-*, 5P(") = 5p(")(t)e"*, (25) 

and 

= zi"\t)e''\ (26) 

where (jy^'^^Zjt), zi"\t) and 6P^°'\t) are periodic functions of time with period 27r/a;, and 
a is the Floquet exponent which determines stability. When 9fJ[cr] becomes positive, pertur- 
bations grow exponentially and the system becomes unstable. The remainder of this paper 
concerns the calculation of the Floquet exponent. 

Upon applying Floquet theory, the resulting equations of motion in terms of the stream 
function 0^"-* and the interface Zg are 

p{dt + a){dl-qiy^^'^^ + ipqJ:\dl - g^)0(") - ipqA^-\dlut^) 

pt (27) 

= / dt' G^''\t - t'){dl - g2)20(")e-(*'-*), 

J — oo 

and 

{dt + a)zs + iqxU^b\z=QZs = -«gx0^"''U=o- (28) 
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In addition, the no-slip and no-penetration boundary conditions at the planes become 

<^(2)(2 = H) = d,^^^\z = H) = ^^^\z = -h) = d,^^^\z = -h) = 0. (29) 
and the interface conditions are 

0(2)(^ = O) = 0W(^ = O), (30) 

and 

(z = 0)z, + d,^^^\z = 0) = d,ul^\z = 0)zs + d,^^'\z = 0). (31) 

For small Zs we have rii = [d^ZsjO, —1) and ti = [1,0, d^^Zs)- Then, the continuity of the 
tangential shear stresses implies that 

<5a(2) 1,^0 = 5a« 1,^0, (32) 

where 

Sa^) = I dt' G^''\t - t'){dl + g2)0H(t')e-(*'-*). (33) 

J — oo 

Note that the pressure term in the stress tensor vanishes because SijUitj = 0, and the 
contribution of the base flow to this boundary condition is not present because the base flow 
is continuous at 2 = 0, and there is no density difference between two fluids. If the densities 
of two fluids are different, one cannot neglect the contribution of base flow. On the other 
hand, the normal component of the stress tensor must be balanced with the surface tension 
as 

where 

-j dt'G("H^-O(5'-3g')a,0(")e'^(*'-*). (35) 

J — oo 

In order to make analytic progress, we focus only on long-wavelength perturbations. For 
small qx, we expand the perturbations equations (25) and (26) in q^ such that 

0(")(^,t) = 4")(^,t) + 0S")(z,t)g. + ^t\z,t)ql + . . . , (36) 



7.=n 



iqlTzs 



(34) 



Zs{t) = Zsfi(t) + Zs,l{t)qx + Zs,2{t)ql + • • • , 



(37) 
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a = (To + (Jiq^; + a2ql + . . . , (38) 

where at each order in the coefficients of (j)^°'\z^ t) and Zs(t) are time-periodic functions of 
2ti /uj, whereas the coefficients of a are independent of time. Next, after substituting these 
perturbations into Eqs. (27) and (28) we solve the set of dynamic equations order by order 
in q^. 

Before delving into the calculation, we can determine ffist from the equation of motion 
of the interface, Eq. (28), at which order in the ffist non- vanishing Floquet exponent 
appears, and its functional form in the base flow uj, and the perturbed stream function. At 
0{1) in q^ the equation of motion of the interface becomes 

aoZsfi + dtZsfi = (39) 

According to Floquet theory Zgfl must be a periodic function of time. Then it follows from 
this equations that (Tq = and Zs.o = Coj constant. 
At 0{q.j:), with ctq = and = (q, we have 

(TiCo + dtZs^i + iUb\z=oCo = -i(po\z=o- (40) 

The time- dependence of 0o"^ (t) can be inferred from the continuity condition of the x com- 
ponent velocity at the interface Eq. (31). Since Zg^o = Co and Ub ~ e^*"^*, 0o is also periodic 
in time, and we take 

4-\z,t)=4-l{z)e-^-' + 4%)e^-K (41) 

Then, since Zg^i must be a periodic function in time because of Floquet theory, it should 
follow again from Eq. (40) that ai = and 

dtZs,i = -i(ub\z=oCo + 0o|^=o), (42) 

from which we find 

e-^^*-c.c., (43) 

2=0 

where c.c. stands for complex conjugation. Note that [^^^i]* = — 2s,i- In calculating Zg^i it 
does not matter whether we use the base flow and stream function of fluid 1 or the fluid 2 
because of the continuity condition at the interface. Since ai = 0, the ffist non-vanishing 
contribution to the Floquet exponent appears at 0{ql). 



bj . 



0,- 
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Consider then the problem at 0{ql), with ao = di = and Zg^ = Co- We have 



?! \z=0- 



(44) 



Similarly to the case at 0{qx) the time-dependence of 0^"^ is determined from the boundary 
condition Eq. (31). In this case, since Zg^i is a periodic function of time rather than a 



constant, and dzU 



(")/ 



0)zs^i has both a non periodic part and a time periodic part, the 



perturbed stream function at 0{qx) is also a sum of non-periodic and periodic parts as 
^["^p. Thus, from Eq. (44) we find the periodic part 



T{a) _ , (a) 

~ V^l.pcriodic 



dtz. 



s,2 



\z=0Zs,l + (pi \z=0 



periodic 



(45) 



and the first non-vanishing contribution to the Floquet exponent is a = (720'^ where 



0-2 



CoL 



u 



NP 



with 



u 



= 0^8,1 



NP 



1 



u 



(1,^(1) 



c.c 



J z=0 



(46) 



(47) 



Therefore, in order to determine the stability of the configuration under study, we need the 



first order distortion of interface Zg^i and the non-periodic part of 



0). 



We show the details of the calculation of the required contributions in Appendix A. A 
complete analytic solution can be found in the limits of small Reynolds number and small 
elasticity. We mention that a complete analytic solution for arbitrary Reynolds number can 
also be found; however, it is too complicated to present here, and the Floquet eigenvalue 
needs to be evaluated numerically. As the limit of relevance in the copolymer case is that 
of vanishing Reynolds number, we do not pursue this general case here. Nevertheless, Ap- 
pendix B describes the steps necessary to find the Floquet eigenvalue for arbitrary Reynolds 
numbers and elasticity. 

For convenience, we define a complex Reynolds number 



(48) 



G*{")(a;) 1 

where ^f^") = G'(°)(a;)/G"(")(a;), and Re^'^^ = pu'^h'^ /G"^°'\uj) are the conventional Reynolds 
numbers of fluids 1 and 2 with Re = Re^^\ In addition, we define n*^") = G"^^\uj) / G"^"\uj) 
with n = n^'^\ The limit of small Reynolds number implies that both Re^^'^ and Re^"^^ H"^ / h"^ 
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are small. Now we solve the equations of motion (27) for the stream function order by order 
in Qx by taking the expansion in Re and g^"'^ as 

= 4t,oi') + Re^^iz) + 0{Re', {g^'^^Y), (49) 

and 

'^Sp(^) = <^Sp,o(^) + ^e0Sp,i(^) + ^(^^'' i9^''^n (50) 

where 0o"-,o('^) ^^"^ '^i'npo(-^) linearly proportional to (^^"^ and 0o"l^i(2;) and 4>^^-^pi{z) 
do not depend on g^"'^ (See Appendix A for details). For long wavelength perturbations it 
is found that the surface tension T does not enter in the calculation because it appears at 
third order in q^. 

From Eq. (46) with Eqs. (A15), (A29), (A62) and (A66), we find for small Re and g^"^ 

(J2 = ^ei + a^isRe + 0{Re\ {g^'^^f), (51) 

where the explicit expressions of Uci and a^\s are given in Eqs. (A72) and (A73). The term dei 
that remains as Re — is a contribution of elastic origin. We find that ad (x {n — l){nH'^ — 
h'^)(g^'^^ — g^^^), where the proportionality factor is positive. This term vanishes if either the 
elasticity (g^"^^ = g^^^) or viscosity (n = 1) of both fluids is the same. In short, we find a long 
wavelength instability due to elasticity stratification when {n — l){nH'^ — h?){g^'^^ ~9^^^) > 0. 
Figure 2 shows a typical stability diagram. 

In addition, there is a contribution purely of viscous origin that appears linearly in Re, 
(Tvis- It agrees with the earlier calculation of Ref. 27 in which the constitutive relation of the 
phases was assumed to be uniaxial but Newtonian. In this latter case, a long wavelength 
interfacial instability occurs whenever the thinner layer is more viscous. 

An important difference between our calculation and the case of Newtonian layers 
is the fact that here the up-down symmetry of the configuration is lost: a^i changes 
sign upon reversing the position of the fluid layers and their thicknesses. The sign of 
a^i{H,h,G<^\uj),G<^\uj)) is opposite to the sign of a,i{h, H,G<^\uj),G<^\uj)). There- 
fore, stability depends on which layer is adjacent to the boundary being sheared. 

The instability mechanism in the Newtonian case arises from the out of phase evolution of 
the vorticity in both fluids'^^''^'^. Small periodic distortions of the interface induce a vorticity 
perturbation of the same sign across the interface, but its sign alternates at each trough 
and peak. With negligible inertia, the vorticities are advected by the shear to create out 
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FIG. 2: Stability Diagram of the two layer configuration. The vertical neutral stability line is 
set by G'^^\uj) = nG'^'^\uj) while the horizontal line h = y/nH. For example, we have taken: a) 
n = 5/6 and b) n = 6/5. 

of phase components of vorticity because of the viscosity contrast. As a result, vorticities 
of opposite sign at adjacent troughs and peaks get close together, and induce a vertical 
motion of the interface. The interfacial instability due to elasticity contrast found here has 
a similar origin in the out of phase vorticity as in the viscous case. The requirement that 
a viscosity contrast must exist {n ^ 1) for instability to occur ensures that there should 
exist an imbalance of vorticity across the interface. The essential difference between the 
viscous and elastic cases is that the out of phase motion of vorticity is driven in our case 
by the elasticity response. As can be seen from the complex wavevector of the base state. 
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Eq. (21), elasticity induces a phase shift relative to the driving shear flow. Similarly, the 
elastic contrast of two fluids generates out of phase components of perturbed vorticity, and 
drives the interfacial instability. 

IV. DISCUSSION 

We discuss in this section the possible implications of our flndings on orientation selection 
of block copolymers under shear. Following a quench of a large aspect ratio conflguration 
(the aspect ratio being the ratio between the lateral dimension of the system and the lamel- 
lar spacing), a spontaneous distribution of locally ordered lamellar domains with random 
orientation results that coarsens with time^''. Under the application of a shear, the wave- 
length of those domains with orientation that has a transverse component (q || v) changes, 
fact that greatly reduces their range of stability'^^. Therefore, after some initial transient, 
it is expected that a majority of domains would be oriented either parallel or perpendicular 
to the shear. Within the order parameter description used here, the local free energy of 
domains along the two orientations is the same, even under shear (since q ■ v = 0). Hence, 
they are in coexistence. 

As noted above, however, viscosity or elasticity of A or B rich regions are typically 
different. In PS-PI the glass transition temperature of the two blocks is quite different, and 
hence their mechanical response within the copolymer is expected to be different. However, 
given the small size of the blocks (on the order of tens of nm) and the large average viscosity of 
polymer melts, flows at the lamellar scale are expected to be negligible under most conditions 
of interest Furthermore, their contribution in slightly distorted lamellar conflgurations is 
of second order in the distortion, and do not contribute appreciably to their relaxation or 
stability^*^. Therefore viscoleastic contrast between blocks seems insufficient to drive short 
scale flows that could account for differences between parallel and perpendicular orientations 
under shear. 

Whereas short scale flows are strongly damped, this is not so for flows that couple to 
long range distortions of a lamellar conflguration^^, or flows at the scale of the characteristic 
domain size in a polycrystalline conflguration. The analysis undertaken here has aimed at 
establishing the relative dynamic stability of parallel and perpendicular domains when long 
wavelength hydrodynamic flows are allowed. We have found that there exists an interfacial 
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instability of hydrodynamic origin in that limit due to viscoelastic contrast between the two 
phases even for zero Reynolds number. 

Domain coarsening and orientation selection are intrinsically nonlinear problems whereas 
we have only addressed here the linear stability of a particular configuration. The two 
can be qualitatively related as follows: Assume that after some transient there will be a 
preponderance of parallel and perpendicular domains in a large aspect ratio sample under 
oscillatory shear. The instability described will set in whenever two such domains meet, and 
the conditions for instability are satisfied (viscoelastic contrast) and low enough wavelength 
(or equivalently, a large enough characteristic domain size). Once any boundary between 
two such domains becomes unstable, a secondary flow is established which is normal to the 
boundary (shown schematically in Fig. 3). Advection of order parameter by this secondary 
flow distorts the parallel domain (v || q) but not the perpendicular (v ± q). The region with 
locally distorted parallel lamellae will have a higher free energy than the adjacent region of 
perpendicular lamellae, and boundary motion will follow to reduce the free energy imbalance. 
As a consequence, we would expect systematic motion of the boundary towards the parallel 
region every time there is a boundary instability. Although we cannot predict the nonlinear 
evolution of the boundary, and hence the evolution of an ensemble of domains in a large 
aspect ratio sample, this instability can provide for a dynamical selection mechanism that 
favors the perpendicular orientation whenever there is viscoelastic contrast between the two 
orientations. This result is generally consistent with the experimental survey of Ref. 19. 

Three main issues remain unresolved. First, the argument presented does not account 
for the preponderance of the parallel orientation when there is no contrast. Second, as 
shown in Fig. 2, there are regions of parameter space in which the parallel-perpendicular 
boundary is stable. Therefore, whether a distribution of domains in a large aspect ratio 
sample would ultimately coarsen or evolve into coexistence cannot be conclusively answered 
by our analysis. Third, although there is ample evidence of viscoelastic contrast between 
differently oriented lamellae, there is no experimental study that we are aware of that has 
directly probed the uniaxial hydrodynamic properties of a bulk lamellar phase. 
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FIG. 3: Orientation selection. The secondary flow induced by the unstable configuration is indi- 
cated by light blue arrows. This flow distorts the parallel lamellae without afl'ecting the perpen- 
dicular lamellae. 
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Appendix A: Expansion results for low Reynolds number and small elasticity 

1. 0(1) in 

At zeroth order in the momentum conservation equation Eq. (27) reduces to, with 
(To = 0, 

pd,dr^t^ = f dt' G^'^\t - t')dt4-l (Al) 

J — oo 

We have shown earher that (f)^^ is time-periodic in 27r/uj, and given by Eq (41). For e~*^* 
mode, we have 



For small Re, with Eq. (49), the above equation (Al) becomes at zeroth order in Re 

d'MlM = 0, (A3) 

and at first order in Re 

1 + 



d'M, = -^T^d!4:L- (A4) 



First, the general solution of Eq. (A3) is 

4%i^) = 4°^ + Blr^z + Ct^z' + D}r^z\ (A5) 

The coefficients are determined by applying the boundary conditions Eqs. (29) - (34). At 
zeroth order in Re the boundary conditions are 



,o(^ = H) = 0, (A6) 



dJ^olA^ = H) = 0, (A7) 

4o,-,oi^ = -h) = 0, (A8) 

d.^olA^ = -h) = 0, (A9) 

0ilo(^ = O) = 0ilo(^ = O), (AlO) 

d,u%iz = 0)zs,o + 9,4'-,o(^ = 0) = d,u%iz = 0)zs,o + 9.4!-,o(^ = 0), (All) 
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G<'\~u:)d'J^l,{z = 0) = G<'\-u;)d'Ml^,{z = 0), (A12) 

G<'\-u:)dl^^^Uz = 0) = G<'\-uj)dl4Uz = 0), (A13) 
where the base flow is also expanded in Re 

u'-^\z) = u^^liz) + u^^\{z)Re + 0{Re\ (g^^^f). (A14) 

In terms of the coefficients of 0q'^1 q this set of boundary condition can be expressed as a 

8x8 matrix. We solve this linear system for the coefficients of 0q"1 q{z) for small elasticity 
gia) ^ obtaining 



.(1) _ .(2) 

/1q — /1q 



{H + h)h^{n-l)nH^CoU 



X 



^0 



{h + nH){'n?m + AnH^h + QnH'^h'^ + AnK^H + /i^) 

-5n'^Hh^ - 2h^n + h^) (A15) 
{AnH^ + h^ + 3nH^h){n - l)hCoU 



X 



Bo 



,(2) _ 



2{h + nH){n'^m + AnH^h + QnH'^h'^ + Anh^H + /i^) 

n[/CoM^7^^^-^^^^) 

\{h + nH)^{n^H^ + AnH^h + 671/72/^2 + 4^/i3^ + /^4)2 

x{-5Hh'^ - 3H'^h^ + 4n^i/S - /i^ - 73n^H^h^ + 2n^H^h^ 

-20nH^h^ - Qn^H'^h - lOH^nh* - bbH^n^h^ - AbH^n^h^ 

+Ah%'^H^ - lAh^nH'^ - 9n^H^h^ + Sn'^H^h - 20/ii7V - 8n^H^) (A16) 

{H^n + Ah^ + 3h^H){n - l)nHCoU 
2{h + nH){n^H^ + AnH-% + Qnmh^ + Anh^H + /i^) 

^^2(/i + nHfin'^H^ + 4nif3/i + 6^1/2/^2 + 4^/i3^ + h'^f ^ 
X (n^i/^ + br/H^'h + ?,n^H^h^ + lAn^H'^h'^ + 20n^H^h^ 

+10n^H'^h'^ - An^H^h^ + 9n^h^h^ + Abn^H^h^ + 73n^H63h^ 

+55n^H^h^ + 20n^Hh'^ - 2nH^h^ + QnHh'^ + 8nh^ - 3HK^ - Ah^) (A17) 
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(1) 



n 



(2) 



{h + nH){n^m + AnH^h + QnE^h"^ + Anh^H + /i^) 



+^ 



(A18) 



(1) 



(2) 



{h^ - nH^){n - l)CoU 



X 



,3 ^4^2 , uh^H^n^ + 2hV*H' + 7n^H'h^ 
-lOnH^h^ - n^H'^ - AH^h^) 



At 0{Re) the general solution to Eq. (A4) is given by 



(a) 



r{°)^2 I r(")-v3 



^0 ^4^:^^5 



(A19) 



(A20) 



where we have used Eq. (A5) for the inhomogeneous term of Eq. (A4). The boundary 
conditions at 0{Re) are 

= H) = 

= H) = 

= -h) = 
-.(1) - ^h) = 



^(2) 



,,{z = 0) = 0« ,(2; = 0) 



dj^\iz = 0)zs,o + = 0) = djl\iz = 0%,o + = 0) 

= 0) = G<'\-u:)d'MUz = 0) 
G*(^)(-a;)af4^1,,(z = 0) = = 0) 



(A21) 
(A22) 
(A23) 
(A24) 
(A25) 
(A26) 
(A27) 
(A28) 
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which can also be written as a 8x8 matrix for the coefficients of (/)q"1 ^. The solution is 

. (n - l)nUCoH'^ 

~ \0{h + nHy{n^H^ + AnH^h + QnH^h^ + Anh^H + h'^f ^ 

x{8n^H^ + 9n^H'^h + QSn^H^h + lAOn^H^h^ + mH^^K^ + lb2n^H^h'^ 

+120n^H'^h'^ + 323n^H^h^ + 370n^H^h'^ + 173n^H^h^ - UOnH'^h^ 

-82h^nH - lASnH'^h^ - QOnH%^ - 22h^ - 21h^H) (A29) 

^(1) ^ n(n - l)HUCo ^ 

° 120/i2(/i + nHY{n'^m + 4ni73/i + QnH'^h'^ + 4^/^3/7 + /i4)2 

x(10n^if^° + SOH'^n^h^ + 93H^n^h^ + 78n^H^h + leOifV/i^ 

+477//'^n^/i^ + 128n^H^h^ + + 305i/V/i^ + 2l2H%^h^ 

+588H^n^h'^ + lOOSi/S^/^^ ^ 817H^n^h^ - 535h^nH^ - A97H^nK^ 
-328h^nH - l8Qnh^H^ - 88h^° - 63h^H) (A30) 

p(2) ^ . - l)^Co 

° 120/i(/i + nHy{nm^ + AnH'^h + QnH^h^ + Anh^H + /i4)2 

X {32n^H^^ + 27H\^h + 272n^H^h + SlSi/^n^/i^ + AQ3H^n^h^ 
+180H^n^h^ + A80n^H^h^ + 1068i/V/i^ + 1123/7V/i^ + 397H^n^h^ 
-28H^n^h^ - 32m^n^h* - bAOH^n^h^ - 543i/V/i^ - UbH^n^h'^ 
+8h^n^H^ - IQOH^nh^ - lA7h^nH^ - 72/i^nif - 20/i^°) (A31) 



.(1) _ }_ j{2) 



{n - l)f/Co 

X 



30/i2(/i + nHf{n^m + 4ni73/i + QnH'^h^ + 4n/i3/J + /i4)2 
x(4n^iJ^° + 20ifV/i^ + 3An^H% + 27 H^n^h^ - 20H%^h^ - IQH^n^h^ 

+37if V/i^ - 8lH^n^h^ + QOn^H^h'^ + 22H^n^h'^ - 80H%'^h^ 

-9QH^n^h^ - 7QH^n^h^ - A8h^nH^ - Alh^nH - 55H^nh^ - (A32) 

° 120/i2(/i + nHf{n^H^ + 4ni73/i + 6^1/2/^2 + ^nh^H + /i4)2 

X (9n^i7^ + ll3n'^hH^ + lOGn^/i^i/^ + QOn'^h^H^ + 32An^h^H'^ 

+b22n^h^H^ + 440n^/i^i/^ + GAn^h^H^ + imn^h^H^ + 176n2/i^i7^ 

+160n^/i^i7^ - A2n^h^H^ + 36n^h'^H^ + lAOnh'^H^ + UAnh'^H^ 

+127nh^H + 31/i^) (A33) 
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(2) ^ „ {n-l)nUCo ^ 

° 120h^{h + nHY{n^H^ + AnH^h + QnE^K^ + 4n/i3i7 + /i4)2 

+U2n^h^H'^ + 20An^h'^H^ - mnh^H^ - Unh'^H^ - 37nh^H - 21h^) (A34) 



2. 



At 0{qx), with ctq = o"i = the momentum conservation equation becomes 

pd.d'J-^ - t dt' G^-\t - t')dt^[-^ = zp4"^(5'4"^) - ^Pu^-^d'A-^ (A35) 



Since we are interested in the time-independent non-periodic part of (j)^^^ in order to ob- 
tain the non- vanishing Floquet exponent a2, dt4>i,Np = and the momentum conservation 
equation reduces to 

(A36) 



P 



NP - 



NP 



where the definition of the steady state viscosity is used 

j'OO 

r^H = / dtG^"'\t), 
Jo 

and the non-periodic part of the RHS is given by 



(A37) 



NP 



. (A38) 



Again, in the same way to calculate the solution at (9(1) in we take a solution expanding 
in small Re as Eq. (50). Thus, at first two orders in Re it is the case that 



d'MU^) = 0, 



47(a) 



n 



(a) 



The general solutions of Eqs. (A39) and (A40) are 



a,NP,o' 



z-V^O,+,0 



(°) -.2 , n(°) -.3 



A 



{a) 

^4 
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^^1 -,5 



120 



A 



2 / 



360 



(A39) 
(A40) 

(A41) 
(A42) 



24 



where 



AS") = ^f[uL")(0)]*cJ")+c.c, 



(A43) 



A 



(a) 



([a.uL"i(0)]*Ci") + 3[nL"i(0)]*D(") + c.c) (A44) 

(A45) 



Similarly to the calculation at 0{1) in q^, the boundary conditions for ^i'npoI-^) 



.(z = if) = 0, (A46) 



^.0Sp,o(^ = ^) = 0, 
'^Sp.oI^ = -h) = 0, 

0Sp,o(^ = 0) = 4>ilp,o{z = 0), 



^^^i^NP o(^ = 0) - "^^^S NP 0(2^ = 0) = - S^n^l^^o - = 0) 2;^,o^L\o + <Pol,o 

— C.C, 



and for ^["np i(^) 



^z</'Sp,o(^ = 0) = na^0i^ipo(2; = 0), 
5'0?ko(^ = O)=^^'<^Sp,o(^ = O)' 

0;'ip,,(2; = iJ)=O, 
^.0Sp.i(^ = ^) = O, 



^l,NP,l' 



-/i) = 0, 



^NP,i(^ = 0) = <^S!ip,i(^ = 0)' 



(A47) 

(A48) 
(A49) 
(A50) 

- 2=0 
(A51) 

(A52) 
(A53) 

(A54) 
(A55) 
(A56) 
(A57) 
(A58) 



dz(pflpAz = 0)-d,4>\"LpJz = 0) = - 



(1) 



1 



9zU% 



z=0 



(1) 

Zs,0U_i 



^0,-,l 



z=0 



+ - 

00 

—c.c. 



J z=0 I 



7 ?/(^) 



J 2=0 

(A59) 
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d'MpA^ = 0) = nd'M^Jz = 0), (A60) 



-c.c, (A61) 
can also be separately written as a 8 x 8 matrix. The solutions are 

A['^ = A? 

_ . h^n^U\oH\H + h){g^^^ - g^''^) 

~ ^u{h + nHf{n^H^ + AnH^h + QnH'^h'^ + Anh^H + h'^f ^ 

x{12nH^h^ + Un'^h^H'^ + Sh^nH + IQnH^h^ + 2n'^H^h^ + ISn'^H^h'^ - 2h^H^ 

- 2H^h^ + n^H^ + 2n^i7^/i + An^H'^h + 2n^H'^h + AQn^H^h^ + ASn'^H^h^ 

+8n^H^h^ + Sn^H^h'^ + IQn^H^h'^) (A62) 

(1) ^ (4g3^ + + 3ng2/i)/^ ^(2) 
^ n{H'M + Ah^ + 3h'^H)H ^ 



x{8H^n^h^ + i/^n^ + TOi/^n^/i^ + 2n^H^h^ + UH^nh^ + ASH^n^h^ 
+18H%^h'^ + 2hH'^n^ + SH'^n^h^ + AGn^H^h^ + AH'^n^h - 2h^H^ 

+ 2n'^H'^h + lQh%H^ - 2h^H^ + lAh%^H^ + Sh^nH) (A63) 



. n[/\o/i(ij'=^n + /i3)((7(2) -^(1)) 

x{8H^n^h^ + H^n^ + 70H^n^h^ + 2n^H^h'^ + l2H^nh^ + ASH^n^h^ 

+ 2r2^i7^/i + lQh%H^ - 2h^H^ + 14/16^2/72 + g/i^ni/) (A64) 
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n 



2uj{h + nHy{n^m + AnH^h + QnH^K^ + 4n/i3/f + h^f 

+AH'^n^h - 2h^H^ + + 2n^H'^h + 16h%H^ - 2h^H^ 

+Uh^n^H'^ + 8/i^niZ) (A65) 



. n{n - l)CoU^H^ 

^ 120uj{h + nHf{n^H^ + AnH^h + GniZ^/iS + 4^/i3^ + /^4)3 ^ 

X {-22Ah^n^H^ + 766h^%H^ + SOAnh^H^ + 388Ah'^n^H^ + 43/1^2 ^ 22Q8h^n^H'^ 

632hVH^ - 32H^h'^n^ - IQH^h^n^ + 7Q0h^H^n - 12h^H%^ + Qh^H^n^ 

-Uh^n^H'^ + 332h^^nH - AOh^n'^H^ + HQh^H^n^ + 2040/i^i7^n2 - Q8h^n^H^ 

+A027h^n^H^ - 288h'n^H^ + An^hH^^ + 22n^h'^H^'^ + 9An^h^H^ + IQOn^h^H^ 

n^H^^ + 38h^^H + 8n^hH^^ + 2An%'^H^^ + hl2n'^h^H^) (A66) 

120/iw(/i + nHY{n?H^ + 4niZ3/i + 6ni72/i2 + 4^/i3^ + /i4)3 
x{20h'^H^n^ + A52H^h^n^ + 29H^'^h%^ + 38AH^'^h^n^ 

+2372H^h^n^ + Qn^hH^^ + 2A2H'^h^^n + 307h^'^nH^ + 5046iZ^/i^n3 

+1488/i^^n2i/3 + lOSOiJ^/i'^n^ + 103i/^°/i^r2^ + 1340/i^°iZV + 788iJ^/iV 

+928H^h^n^ + 58iJ^/i V + A8n^h^H^^ + IGGi/^^/i^n^ - 20n^h^H^^ 

+16n'^hH^^ + IbAh^^nH + 8820H^h^n^ + 292h^\^H^ + 2A0h^H%'^ 

+2n'^H^^ + 2789H^hV - 2Ah*HV + lUOH^h^n^ - 48/i^iJV 

+20/i^^ + 828H^h^n^ + AA32H^h^n^ + 8536if^/i^n^) (A67) 
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. n{n - l)CoU^H 

" ^ 120hu{h + nHf{n'^H^ + AnH^h + QnH'^h'^ + Anh^H + h'^f ^ 

-QQAnh^'^H + lOn'^hH^'^ - 1337nh^°H^ - 1020n^h^H^ - 2A00n'^h^H'^ 
-355Qn'^h'^H^ - 2782n^h'^H'^ + Un^hH^^ - UQAn'^h^^H'^ + IQn^h^H^ 
-nOOn^h^H^ - A29n^h^H^ - 7Q8n^h^H^ - 3b2n^h^H^ - 77A2n^h^H^ 
-822n%^H^ - 15A0n^h^H^ + 23n^h^H^^ - UOQnh^^H^ - A5Qnh^H^ 
-USAn^h^^H^ - mn%^H^^ - 86/i^^ + 32n%^H^^ - 34n^h^H^° 
-32n^h^H^ + 2n^H^^ - 57h^^H) (A68) 

n 

. {n - l)CoU' ^ 

^ 120uh'^{h + nHy{n^H^ + AnH^h + QnH'^h^ + Anh^H + h^^ ^ 
X {-8H^^'h^n^ + n'H^^ + eOOSif^/i^n^ - An^h^H^^ + 2464if^/i^n^ 

+157QH^h^n^ + 270AH^h%^ + A80h^H^n^ + 6093H'^h^n^ + AOOH^h^^n 

+357AHVn^ + 2An^h'^H^'^ + l22AH^h^n'^ + 832h^^H%^ + 3Q8m^h^n^ 

+8n^hH^^ + 52Ah^%H'^ + 2148/i"n2i73 + 3S2h^'^nH + 336H^°h^n^ 

+5QH^^h^n^ + 1040if^/i*^n2 + 800H^h^n^ + Q32hVH'^ + 308H^h^n^ 

+3552H^h^n^ + + 43/i^^ + 2391i/*^/i^n^ + 182i7i°/i^n^) (A69) 

x{2n'^H^^ - 252nh^^H^ - 29h^^ - IQn^h^H^ - lOAGn^h^H^ 
-33Anh^^H^ - 90n^h^H^^ - 68An^h^H^ - A36n^h^^H'^ - An^h^H^^ 
-llQAn^h^H^ - l89Qn^h'H^ + 15n^hH^^ - 18Qn^h^H^ - 900n'^h'^H^ 
-2060n^h^H'^ - 1296n^h^0H^ - lOOn^h^H^ - lA92n^h^H^ - 22Qnh^^H 
-Ib&n^h^H* - 176n^h^H^^ - 2059n^h^H^ - 57Qn^h^H^ + 8n^h^H^^ 
-320n^h^H^ - 8n^h^H^^ - 13Q7n^h^H^ - 2A0n'^h^H^ - lAOAn^h^H'^) (A70) 
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-A52n^h^H^ - 106nh^^H + 2804^^/1^/7^ + 144n^/i^/7^° - 389n^h^H^ 
-256n^h'^H^ + 108n^/i2i/^^ + Un^h^H^^ + 1583n^/i^i/^ + lOAOn^h^H^ 
+223Qn^h^H'^) (ATI) 

We now obtain the Floquet exponent Eq. (51) using Eqs. (A15), (A29), (A62) and (A66), 
resulting in 

_ U^h^H^H + h){n- l){nH^ - h^){nH^ + h^){g^^'> - g^^')) 
~ co{h + nH)'' 



n 



^H^ + 2hn'^H^ + 2H^n^h^ + 2nH^h + AnH^K^ + Ah^nH + 



^ {h^ + GnH^h^ + n'^H^ + Ah^nH + AnH^hf ^^^'^^ 



and 



U'^H'^n{n - l){Ahn^H^ + n^H^ + + Ah^H + QnH'^h^) 
^'''^ ~ 120w(/i + nHy{h^ + GnH^h"^ + ^2/74 + 4/^3^^ + AnH^hf ^ 

X (n^i/s + 471=^/7^/1 - 2n^H^h'^ - Ah^n^H^ + Sn'^H'^h'^ + Sn^H^h^ 

-Sh^n^H'' - Sh^^H^ + Ah^nH^ + 2h%H^ - Ah^nH - h^) (A73) 

Appendix B: Expansion for general Reynolds number 

Now we present a general framework to obtain a solution of the equations of motion for 
an arbitrary Reynolds number, but still in the limit of small qx- 

1. 0{1) 

For arbitrary Reynolds number we can find a general solution to Eq. (Al) by taking 
Eq. (41) for ^o""*. For e"*'^* mode the general solution to Eq. (Al) is 

4") (^) = + 5$")^ + Ct^e'""'' + Di;'^e-''^'' (Bl) 
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where fc^"^ is given by Eq. (21). Again, the coefficients of (^^^^{z) are determined by applying 
the boundary conditions, Eqs. (29) - (34) which are 



-kWh 








kWh 



1 



(2) 



5, 



(2) 



(2) 



5, 



(1) 



(1) 



(1) 



(B2) 

(B3) 
(B4) 
(B5) 
(B6) 



(^)^~-0)Co + 5(^)+A:«C'(^)-A;«I)« 

(B7) 
(B8) 

(B9) 



Therefore we obtain 0q") by solving this 8x8 linear system. 



2. Oiq,) 

At 0{qx), we have to solve Eq. (A36) for ^i^^p. If Eq. (Bl) is used in the RHS of 
Eq. (A36), we get 



/ '^z'f'l.N 



NP ~ "^P 



(A;(°))2 _ (A;(")*)^ 



(a) 



-ip(fc(")*)2MP(ii") + 5^")^) - C.C. 



(BIO) 



The solution to the above differential equation has a homogeneous solution ^["np h ^ 



particular solution 0["np,p- First, the particular solution can be found by using the fact that 



dtf{z) = {A + Bz)e 



±kz 



(Bll) 



so that 



and, similarly. 



dt3(z) = Ce 



±kz 



(B12) 



(B13) 
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C 

The resulting particular solution for fluid 1 is 



±kz 



^(1) 
l,NP,p 



u 



A, 



(1) 



5, 



(1) 



+ 



5, 



(1) 



{kW*y 



(1) 



4 



5, 



(1) 



+ 



5, 



(1) 



(A;(i)*)4 (A;(i)*)5 {kW*) 



^kW*{z+h) 



-kW*{z+h) 



4 7^(1) 



I p 



(A;(l))2_(^.(l)*) 



(1) 



(fc(l)*+fc(l))2 



+ 



(1) 



477(1) 



*n2 



X 



(1) 



K* 

(fc(l)*_fc(l))2 



+ 



(1) 



_(fc(i).+fc(i))2 



-c.c. 



(B14) 



(B15) 



and for fluid 2 



l,NP,p 



(2) 



(2) 



5, 



(2) 



(A;(2)*)5 (A;(2)*)^ 



(2) 



+ 4 



(2) 



+ 



5, 



(2) 



(A;(2)*)4 (fc(2)*)5 (A;(2)*) 



4 7^(2) 



i^* G*(2)(^)^(2)* 



A, 



(2) 



(A;(2) 



4 ^° 



(2) 



{m*f (fc(2) 



A, 



(2) 



+ 4 



5, 



(2) 



+ 



5, 



(2) 



(A;(2)*)4 (A;(2)*)5 (A;(2)*) 



+ 



4 7^(2) 



(A;(2))2 _ (A;(2)*)2 
(2). 



+ 



sinh(A;(^)*/i) 

(fc(2)*_fc(2))2 



4r/(2) L 



I ^^(2)* + ^(2))4 ' (^(2)* _ ^(2))4 



(1) 



-(fc(2)*+fc(2))-. 



+ 



(A;(2))2 _ (A;(2)*)2 
^(fc(2).+fc(2))2 



(2) 



g_(fc(2)._fc(2))2 



I ^^(2)* _ ^(2))4 ' (^(2)* +^(2))4 

cosh(A;(^)*/i) 



(A;(2)* + A;(2))4 (^^(2)* _ ^(2))^ 



,(fc(2)*+fc(2))2 



.(fc(2)*_fc(2))2 



-C.C. 



(1) 

I ^^(2)* _ ^(2))4 (^,(2)* + ^(2))4 

(B16) 
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Note that 4>^^^p p{z) is purely imaginary. On the other hand, the homogeneous solution is 



(B17) 



The eight coefficients of ^["np h determined by imposing the eight boundary conditions. 



^(2) 



d.^?i,P,Hi^ = H) + d.4>\%p/z = H) = 



z = H) + ^pf^pjz = H) = 



(2) 



^Sp,.(^ = 0) + 0Sp,,(^ = 0) = 0« p,,(z = 0) + 0;;ip,^(z = 0) 



dzu''^\z = 0)zs,i 



9. 



NP 



/,(2) 

a,NP,/i 



c(2) 
l,NP,p 



+ dz 



c(2) 

a,NP,/i 



+ 



l,NP,p 



z=0 



NP 



2=0 

r'l.NP.h, 



+ 



^(1) 
l,NP,p 



2 = 



^(1) 

a,NP,/i 



+ 



ip 



J NP,2=0 



7] 



.(2) 

''l,NP,/i 



l,NP,p 
'f'l.NP.p 



2 = 



J 2=0 



ip 



NP,2=0 



V 



(1)9? 



^(1) 

a,NP,h 



l,NP,p 



2 = 



(B18) 

(B19) 
(B20) 
(B21) 
(B22) 

(B23) 
(B24) 
(B25) 



Again, the general solution can be found by solving this 8x8 matrix for the coefficients of 



^l,NP,h. 



In principle, we can get a2 of Eq. (46) from the derived solutions 0i"np,/i ^ij^p p 

-,2 



.r 

so 



u 



{l)l 

b \z- 



NP 



^ ^(") I 



^(") I 
''l,NP,pU=0 



(B26) 



where 



and 



u 



(")| 

6 



=0^s,l 



NP W 



(1)1*^(1) 



0, 







c.c. 



2=0 



(1) 



(B27) 



(B28) 
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lH I 



^ P nJi)*\2 U ^kW*h 
4 7]^^') K* 



A, 



(1) 



-4- 



(1) 



u 



4?7« 
4r/« 



(1) 



, . ^0 



(1) 



*^4 



+ 4 



(^(l)*)£ 



*\2 



4r/W 
-c.c. 



(A;a))2_(fc(i)*)^ 



1 t/ 



(1) 



+ 



(1) 



+ fc(l))4 



(1) 



(1) 



(B29) 



In general, (T2 is too a complicated function of fluid parameters to present explicitly, but 
we can readily evaluate it numerically. However, one should be careful in the numerical 
evaluation as we have observed numerical instabilities in certain limits. For example, in the 
limit of small Reynolds number 3fJ[A;*^"''] becomes too small, and causes singular behavior, 
leaving the numerical result unreliable. In this regard, we decided to work directly on the 
Umit of small Reynolds number as done in Appendix A 



